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[57] ABSTRACT 

A method and a test instrument are provided for measuring 
program clock reference (PCR) jitter in a data stream, such 
as an MPEG-2 transport stream. The test instrument includes 
a line interface circuit for receiving the data stream, a digital 
computer for processing the data stream and a local clock for 
generating timestamps. The digital computer extracts from 
the data stream a PCR sample value which represents a 
count of a system time clock used to generate the data 
stream. The timestamp is obtained from the local clock and 
is associated with the PCR sample value. The digital com- 
puter determines a reconstructed system time clock from a 
set of values including the PCR sample value and the 
timestamp, and previously-received PCR sample values and 
associated timestamps within a window representative of a 
predetermined time interval. The window is shifted as new 
PCR sample values are received. The digital computer 
determines a PCR jitter value for the PCR sample value by 
determining a deviation of the PCR sample value from the 
reconstructed system time clock. The PCR jitter value may 
be displayed on a display unit, preferably in a histogram. 

20 Claims, 8 Drawing Sheets 
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METHOD AND APPARATUS FOR PCR data stream containing digitally encoded information and 

JITTER MEASUREMENT IN AN MPEG-2 PCR samples are provided. The data stream typically com- 

TRANSPORT STREAM USING SLIDING prises an MPEG-2 transport stream. The method of the 

WINDOW invention comprises the steps of providing a test instrument 

CROSS REFERENCE TO RELATED 5 ^chiding a line interface circuit for receiving a data stream, 

APPLICATION a ^S^ 3 ' computer connected to the line interface circuit for 

processing the data stream, and a local clock for generating 

This application is a continuation-in-part of application timestamps. The digital computer extracts from the data 

Ser. No. 08/632,036, filed Apr. 12, 1996, now U.S. Pat. No. stream a PCR sample value which represents a count of a 

5,774,497. 3Q system time clock used to generate the data stream. A 

ptPinnPTHP iiMVPNTtnM timestamp, which represents a time of arrival of the PCR 

ntLUUr iHciwviiwiiuiN sample va]uej k obtained from the bca] clock and is 

This invention relates to analysis of data streams trans- associated with the PCR sample value. The digital computer 

mitted through a data communication channel and, more determines a reconstructed system time clock from a set of 

particularly, to test instruments and methods for measuring values comprising the PCR sample value and the timestamp, 

program clock reference (PCR) jitter in an MPEG-2 trans- 15 and previously-received PCR sample values and associated 

port stream. timestamps within a window representative of a predeter- 

DA^rnniTKm i-nmrrj un/PKmnM mined time interval that terminates with the PCR sample 

BACKGROUND OF THE INVENTION yalue and tfae timestamp , llie digital computer determine ; a 

In the field of data communications, one highly demand- PCR jitter value for the PCR sample value by determining 

ing application is the transmission of full motion video, such 20 a deviation of the PCR sample value from the reconstructed 

as may be required for video-on-demand applications. Vari- system time clock. The PCR jitter value represents a varia- 

ous protocols which involve compression of digitized video tion in the arrival time of the PCR sample value with respect 

information have been proposed. One such protocol is the to the reconstructed system time clock. 

MPEG-2 transport stream, as defined in standard ISO/I EC The reconstructed system time clock is preferably deter- 

13818-1: 1996 promulgated by the ISO/IEC The MPEG-2 25 mined by determining a line of best fit to the set of values 

transport stream involves transmission of video and audio within the window. The slope and the ousel of the line of 

information in transport stream packets of 188 bytes in best fit are calculated. The line of best fit represents the 

length. Each packet includes a header, which contains con- reconstructed time clock. 

trol information, and a payload, which contains video or 'ITie above procedure is preferably repeated for a plurality 

audio information. The MPEG-2 transport stream may carry 30 of PCR sam P le values extracted from the data stream to 

multiple different programs simultaneously. Each packet in P rovide a P Iur ality of PCR J iUer value s. ™* PCR jitter 

the transport stream is associated with a program by a packet values are determined in real time as the PCR sample values 

identifier (PID) contained in the header. The header is of are received, without waiting for all PCR sample values to 

variable length, depending on whether it contains an adap- be received - ™e window is shifted relative to the set of PCR 

tation field. The adaptation field contains control informa- 35 sam P !e values ™ d timestamps for each new PCR sample 

tion that is not necessarily present in every transport stream value rece,v ed, 50 that a predetermined number of the most 

packet recently received sample values is used to determine the line 

One' of the fields that may be contained in the adaptation °.f t best f ^ a slidin S window 15 t0 determine PCR 

field is a program clock reference (PCR). The PCR is J values. 

representation of a system time clock that was used at the « ^ ste P of determining the line of best fit preferably 

source to encode the data in the transport stream packet. The mcludes a recursive Unear regression calculation on the set 

system lime clock typically has a frequency of 27 MHz. °£ PCR sam P le values and timestamps within the window. 

Each PCR field contains a PCR sample value which repre- ^ vanance of the PCR jitter values may also be deter- 

sents a count of the system time clock at the time when the mmed usin S a recursive calculation, 

transport stream packet was encoded for transmission. PCR 45 The invention provides an accurate measurement of PCR 

sample values for a particular program are received with j itter for each PCR sample value received. The PCR jitter 

different transport stream packets having the same PID. A measurement is available in real time as the PCR samples 

series of PCR samples can be used to reconstruct the system are received. The accuracy of the PCR jitter measurement is 

time clock at the decoder end of the communication channel. independent of the number of PCR sample values received. 

The MPEG-2 transport stream may, for example, be 50 Furthermore the PCR jitter measurement quickly adjusts to 

transmitted through asynchronous transfer mode (ATM) a chan S e 10 Parameters, such as a change in network delay 

networks developed for high speed, packet ized digital trans- or c * ock s * ew " 

mission of data, audio and video. One of the parameters that BRIEF DESCRIPTION OF THE DRAWINGS 

is used in characterizing an MPEG-2 transport stream is For a better unders t an ding of the present invention, ref- 

PCR jitter. In an ideal cx)mmunxation channel, all transport 55 . g ^ to ^ ^ d ^ ^ 

stream packets and therefore all PCR samples, are received ^0^^ herein by reference and in which: 

a fixed time after transmission. However, m a real commu- . . . . . ,. r . - t 4 

. , • i * j * \ . t , . . FIG. 1 is a block diagram of an example of a test 
nication channel, variable delays may be introduced by . , - ( . . r . , „™ .f, 
,-rr . . . . t r I ,-rr . . i instrument suitable for implementing PCR jitter measure- 
different channel elements. For example, different transport , ... v . .. & J 

, t c „ . ' i ■ . An menl in accordance with the invention; 

stream packets may follow different network paths in reach- 60 ^„ „ . a „ . .. 

ing the final destination. Variations in arrival times produce c ™ „ 2 ? a flow dia 8 ram of a firsl em bodiment of a routine 

PCR jitter. PCR jitter is described in standard ISO/IEC for PCR J ,tlcr measurement; 

1381 8-9: 1996. 3 is a graph of normalized PCR values as a function 

of normalized PCR arrival times; 

SUMMARY OF THE INVENTION 65 FIG. 4 is a block diagram that illustrates hardware and 

According to the present invention, methods and appara- software components involved in an implementation of the 

tus for measuring program clock reference (PCR) jitter in a invention; 
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FIG. 5 is an example of a histogram showing the fre- 
quency of different PCR jitter values; 

FIG. 6 is a graph of normalized PCR values as a function 
of normalized PCR arrival times, illustrating the change in 
slope of the line of best fit that may result from a change in 
network delay; 

FIG. 7 is a graph of normalized PCR values as a function 
of normalized PCR arrival times, illustrating different lines 
of best fit at different window times; and 

FIG. 8 is a flow diagram of a second embodiment of a 
routine for PCR jitter measurement. 

DETAILED DESCRIPTION 

A block diagram of a test instrument for analyzing a 
MPEG-2 transport stream is shown in FIG. 1. Line interface 
hardware 10 receives an MPEG-2 transport stream, which 
may be transmitted from a remote location. The MPEG-2 
transport stream is typically received on ATM and possibly 
a physical layer, such as SONET. The MPEG-2 transport 
stream data is supplied to a protocol processor 12. The 
protocol processor 12 receives each packet in the MPEG-2 
transport stream and identifies information in the header in 
accordance with the MPEG-2 protocol. The protocol pro- 
cessor 12 also includes a decoder for decoding the payload 
information. Selected information can be extracted from the 
MPEG-2 transport stream by the protocol processor 12. 
Suitable line interface hardware and protocol processors are 
known in the art and are commercially available. For 
example, the line interface hardware may be a Model 
El 697 A, and the protocol processor may be a Model 
E4209B, both sold by Hewlett-Packard Company. 

Selected information from the protocol processor 12 is 
supplied to a host controller 14, which analyzes the selected 
information and provides results to a display unit 20. The 
display unit 20 may, for example, be a video monitor. The 
host controller 14 may be a general purpose computer, such 
as, for example, a type E1499, and may utilize the Unix 
operating system. 

The MPEG-2 transport stream analyzer shown in FIG. 1 
may perform functions such as real time analysis of program 
clock reference statistics, packet identifier statistics, pro- 
gram specific information table monitoring, transport stream 
and global error indicators, header statistics and header field 
matching. The test instrument may also provide triggering 
and pattern matching functions, and decoding and display of 
the video contained in the transport stream. 

As indicated above, the MPEG-2 transport stream is 188 
bytes long and includes a header and a payload. The 
following header information is relevant to the present 
invention. The header includes a two-bit adaptation control 
field, which indicates whether the transport stream packet 
includes an adaptation field. The header also includes a 
13-bil packet identifier (PID) field, which identifies the 
program with which the packet is associated. The MPEG-2 
transport stream may carry multiple programs 
simultaneously, with each program identified by the PID. 
The adaptation field includes a one-bit PCR flag, which 
indicates whether the adaptation field contains a PCR field. 
The program clock reference (PCR) field is a 42-bit field 
coded in two parts, a 33-bit program clock reference base 
and a 9 -bit program clock reference extension. The PCR 
base and the PCR extension are defined by the MPEG-2 
protocol. The PCR indicates the intended time of arrival of 
the byte containing the last bit of the PCR base at the input 
of the system target decoder. The PCR represents the time at 
which the packet was transmitted in units of a system time 
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clock count. PCR samples are usually not included in every 
transport stream packet. 

The present invention provides a method and apparatus 
for measuring PCR jitter in the MPEG-2 transport stream. 
The PCR jitter is defined as the difference between the actual 
arrival time of a PCR sample value and the expected arrival 
time based on system time clock (STC) at the source of the 
transport stream. Before the PCR jitter can be estimated, the 
system time clock must be recovered from the PCR sample 
values in the transport stream. The PCR jitter can then be 
measured as deviations from the system time clock. The 
system time clock is subject to frequency variations, but 
such frequency variations can be ignored for short segments 
of the transport stream. In determining PCR jitter, the system 
time clock frequency is assumed to be constant. 

A flow chart of a first embodiment of a routine for 
determining PCR jitter is shown in FIG. 2. PCR jitter is 
determined by a software routine in the host controller 14 
without requiring special dedicated hardware. The PCR 
jitter measurement routine may, for example, be imple- 
mented using the C programming language. In step 40, a 
PCR sample value is extracted from the transport stream. As 
indicated above, several programs may be carried by the 
MPEG-2 transport stream simultaneously. One or more 
programs may be selected for analysis of PCR jitter. These 
programs are identified by the program identifier, PID, in the 
header of the transport stream packet. Although the instru- 
ment may analyze more than one program simultaneously, 
the analysis of one program is described. The routine is 
repeated for each program of interest. When a transport 
stream packet with the selected PID and containing a PCR 
sample value is received, the PCR sample value is extracted 
from the transport stream and is assigned a timestamp from 
a local clock. The timestamp establishes the arrival time of 
the PCR sample value. The received PCR sample value and 
the timestamp from the local clock are normalized in step 
42. The received PCR sample value is normalized as fol- 
lows: 



PCR-^PCR^ „x3Q0)+PCR^ t J/27 



(1) 



where PCR tew represents the 33-bit PCR base field, PCR^ 
represents the 9-bit PCR extension field and PCR„ repre- 
sents a normalized PCR sample value, given in units of 
microseconds. The timestamp ticks provided by the local 
clock are counts of tens of nanoseconds. Normalized times- 
tamp values TL M are obtained by dividing the raw local clock 
count by 100. Thus, the data associated with each PCR 
sample value includes a normalized PCR sample value, 
PCR M , and an associated normalized timestamp, TL„, both in 
microseconds. 

In order to determine PCR jitter, the system time clock is 
recovered from the PCR sample values, and the PCR jitter 
for each PCR sample value is measured as a deviation of the 
PCR sample value from the recovered system time clock. 
The system time clock may be recovered using a linear 
regression model based on the local timestamps of the 
transport stream packets containing the PCR sample values. 
The timestamp is based on the arrival time of the last byte 
in the transport stream packet containing the PCR. 

After two PCR sample values are received, a line con- 
necting those sample values may be drawn. The line repre- 
sents the first estimation of the system time clock. As further 
PCR samples are received, the least squares linear regression 
model may be used to make successive estimations of the 
line of best fit, which represents the reconstructed system 
time clock. A graph of normalized PCR sample value, PCR„, 
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as a function of normalized timestamp, TL n , should 
approach unity slope when there is no frequency variation in 
the encoder system time clock or the local clock. 

A plot of PCR in microseconds as a function of TL in 
microseconds after several PCR sample values is shown in 
FIG. 3. Normalized PCR samples values, PCR„ and corre- 
sponding normalized timestamps, TL rt , are plotted as points 
60, 62, 64, etc. A line of best fit 70 can be determined for any 
set of two or more sample values. The line of best fit 70 
represents the recovered system time clock from the source 
of the transport stream. The line of best fit 70 is a graph of 
the time at which the packet was transmitted (normalized 
PCR sample value) as a function of PCR arrival time 
(normalized timestamp). The line of best fit 70 should have 
a slope of one or approximately one. The PCR jitter, PJ„ of 
each PCR sample value is the horizontal distance parallel to 
the TL axis between each sample point and the line of best 
fit 70. The PCR jitter represents a time deviation of the PCR 
sample value from the recovered system time clock. 

In order to determine PCR jitter, the line of best fit 70 may 
be first determined using a linear regression model. The line 
of best fit is given by: 

PCR-m„xTL+c„ (2) 

where m„ represents the slope of the line of best fit, and c„ 
represents the offset (the PCR axis intercept) of the line of 
best fit after n samples. The values of slope m„ and offset c„ 
are determined in accordance with the linear regression 
model by the following equations: 



i I TLiPCR,- 
i-O 



(" \( " \ 

t, - I 71, I I X PCRi I 



<x X TL? 
/-O 



"I * n ) 

V - 0 / 



X PCRi-m X 
i=0 i=0 



71, 



On each successive sample (TL„, PCR„), the slope m„ and 
the offset c„ of the fine of best fit are recalculated in 
accordance with Equations (3) and (4). However, it will be 
understood that it is not necessary to repeat each calculation 
in its entirety. Instead, the sums in the expressions for slope 
and offset are saved after each calculation. When each new 
sample (TLj,, PCRJ, is received, the new values are added 
to the appropriate sums in the expressions for slope and 
offset to provide new values of slope and offset. Thus, a new 
line of best fit is determined for each new PCR sample value 
with a minimum of calculation. As discussed below, by 
calculating the line of best fit in real time, the PCR jitter and 
PCR jitter variance can also be calculated in real time, rather 
than waiting until all PCR sample values are received. If it 
is assumed that there is little or no slew in the system time 
clock and the local timestamping clock, the successive 
estimations of slope and offset become stable after an initial 
training period and the mean jitter approaches zero for a 
large number of PCR samples. 

Referring again to FIG. 2, when a sufficient number of 
normalized PCR M values have been received, as determined 
in step 46, the slope m„ and offset c„ of the line of best fit 
are calculated in step 50, as described above. When a 
sufficient number of PCR„ values has not been received, the 
process simply wails for more PCR sample values to be 
received in the MPEG-2 transport stream. Preferably, the 
process waits until the slope and offset are stable within 



defined limits before using those values in the jitter calcu- 
lation. In an alternative approach, a line of best fit can be 
determined after two PCR sample values have been 
received. The line of best fit is then updated for each new 

5 PCR sample value received until a sufficient number of PCR 
sample values have been received to provide a desired 
accuracy in the PCR jitter measurement. In step 52, the PCR 
jitter, PJ„, is determined from the normalized PCR sample 
value, PCR„, normalized timestamp, TL„, and from the 

io slope m n and offset c„ of the line of best fit in accordance 
with the following equation: 
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(5) 



The PCR jitter variance, VPJ„ may optionally be determined 
in step 54 in accordance with the following equation: 



VPJ n - X 



WW 



20 



The values of PCR jitter and PCR jitter variance are mea- 
surements of the physical timing characteristics of a selected 
program in the MPEG-2 transport stream. The PCR jitter 
measurement provides an individual PCR jitter value for 
each PCR sample value received and is available in real time 
shortly after the PCR sample value is received. The PCR 
jitter variance is a statistical representation of the PCR jitter 
values and is recomputed for each new PCR sample value. 
It will be understood that the calculation of PCR jitter 
variance is optional. The measured values of PCR jitter and 
PCR jitter variance may be displayed and/or stored in any 
desired manner. For example, the measured values may be 
displayed on the video display unit 20, may be printed 
and/or may be stored for subsequent analysis. In a preferred 
embodiment, PCR jitter values are presented on the video 
display unit 20 in a histogram, as shown in FIG. 5. In the 
histogram, the horizontal axis is the PCR jitter, PJ„, and the 
vertical axis is frequency of occurrence. Routines for gen- 
erating a histogram display using a set of input values are 
well known to those skilled in the art. As discussed above, 
PCR jitter values are determined in real time. The histogram 
may be updated for each new PCR jitter value determined. 
It will be understood that different display formats can be 
utilized for presenting the PCR jitter measurements. 

A preferred implementation of the present invention is 
described with reference to FIG. 4, which shows portions of 
the protocol processor 12 (FIG. 1) and the host controller 14. 
The host controller 14 provides a PID list to the protocol 
processor 12. The PID list specifies one or more programs in 
the MPEG-2 transport stream that are to be analyzed. As 
indicated above, the packet identifier (PID) associates each 
transport stream packet with a program. The MPEG-2 
transport stream is input to a capture RAM 102 in the 
protocol processor 12 through a PID filter 106. The PID filter 
106 selects transport stream packets specified by the PID list 
and inputs filtered transport stream packets to capture RAM 
102. A local clock 104 provides timestamps that are asso- 
ciated with the transport stream packets stored in the capture 
RAM 102. A microprocessor 108 provides raw PCR and 
timestamp values from capture RAM 102 to a shared RAM 
110. The raw PCR data is accessed by the host controller 14 
for computation of the PCR jitter and PCR jitter variance 
values by software routines in block 112. The computations 
are performed as described above in connection with Equa- 
tions (1) to (6). The measured values may be provided to 
computer display 20 as a histogram and are preferably stored 
by the host controller 14 for subsequent analysis. 
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A histogram may be generated for each program, or PID, Id addition, the sensitivity of the PCR jitter measurement 

being monitored. Each bin in the histogram is defined by a to changes, such as changes in network delay and/or clock 

range of PCR jitter. The user may specify the histogram slew, is improved. A graph of normalized PCR values as a 

range. The bin boundaries are calculated from the range. function of normalized PCR arrival times is shown in FIG. 

Respective bin values are incremented for each matching 5 7. The sliding window technique is illustrated as a window 

value of PCR jitter, PJ„. The PCR jitter variance, VPJ n , and 230 at time tl, a window 232 at time t2 and a window 234 

the number PCR samples may be shown as numerical values at time t3. It is assumed that a change in network delay 

on the display. occurs during time tl. Lines of best fit 240, 242 and 244 are 

The technique for determining PCR jitter described above calculated at times tl, t2 and t3, respectively. As illustrated, 

is generally satisfactory, but has certain drawbacks. In the line of best fit changes in response to the change in 

particular, the summation of product terms and product of network delay. The PCR jitter measurement using the sliding 

summation terms in Equations (3) and (4) are cumulative window approach converges to the ideal line of best fit 

over all PCR sample values received. Since the values of faster, since the ratio of the number of sample values 

these terms increase over time, precision may be lost due to affected by the change in network delay to the number of 

floating point limitations of the computer. Specifically, the sample values prior to the change increases much more 

computer has a fixed number of bits, such as 32 bits, for 15 quickly. As the window slides, it reaches a point where the 

representing the summation terms. As the summation terms line of best fit is determined solely by the sample values after 

increase, the least significant bits are dropped, and precision the change in network delay. 

is lost. The loss of precision reduces the accuracy of the The sliding window approach to PCR jitter measurement 

calculated slope and offset of the line of best fit, which in may be understood by considering sample values (TLj, 

turn compromises the accuracy of the PCR jitter measure- 20 PCRj), (TL^, PCR 2 ), . . . (TL„, PCRJ, (TL^, PCR W+1 ). For 

ment. a window having a size n, the slope and offset of the line of 

Calculating the slope and offset of the line of best fit over best fit are first calculated using the sample values 1 to n. 

all PCR sample values received may compromise the sen- When the sample value n+1 is received, the sample value 

sitivity of the PCR jitter measurement to changes in param- (TLi. PCRJ is discarded, and the slope and offset of the line 

eters associated with the PCR jitter measurement. Examples 25 of test fit are calculated using sample values 2 to n+1. This 

include changes in network delay and/or system clock slew. process is repeated for each new sample value received, so 

The reduction in sensitivity is greatest when such changes that the w ">dow slides with respect to the set of received 

occur after a large number of PCR sample values has been sam P le val f s * J 1 * v™ 6 ™ slid t es ? n f sam P le L value at a 

received by the test instrument. time, such that the new sample value falls within the window 

The situation is illustrated in FIG. 6, which is a graph of 30 and ' he ol u de ^ sam P 1 f e value ls Raided, 

normalized PCR values (PCR) as a function of normalized 11 ™ u ? h ^ °[ ™ ,ve formulas > th K e ™P**">nal 

PCR arrival times (TL). A line of best fit 200 is determined <»>Pk»ty and overhead associated with me calculation of 

using Equations (3) and (4) for a large number of received ! he slo P e and oSsti f th ? lu * 0 1}™} fit > and al f° * e PCR 

PCR sample values. Assume that a change in network delay ^ tter vanance ' ma y be reduced ; ? 1S 15 ach * V( ; d b * d ™S 

occurs following PCR sample value 202. Since the param- 35 »he mathematical equations which express the latest value of 

eters of the line of best fit 200 have been calculated over a CJ f of ^ uanUll f " a of < he u latest sam ^ 

large number of PCR sample values, the sensitivity to the value and lhe P rcviousl y calculated value of the correspond- 

change in network delay is reduced. Due to the change in mg quantity. 

nelwork delay, a belter line of best fit would be the dashed u Smce ' hese ^ ll0 ^ a u re ^by. * * ™ nVen ' e °' L *° 

line 206 shown in FIG. 6. However, due to the weighting of 40 c *™& ^ notation used above. Henceforth, TL, and PCR,. 

the previous PCR sample values, the calculated line of best sha11 ^ «P™*«««« by x, and y,, respectively. Ako, PCR 

fit is likely to be very close to line of best fit 200. As a result, Ji"" shall be represented by j. Using this notation Equation 

the PCR jitter calculation would be less accurate following < 3 > for c ^\ a ^ «»* *V* of the line of best lit can be 

the change in network delay. rewnt,en as follows - 

In order to avoid these drawbacks, a sliding window may 45 / \ / \ m 

be used for calculating the slope and offset of the line of best " / " \( " \ 

fit. Instead of using all the available PCR sample values in 4 " ( Ci x<yi ~ I * f 1 1 ,-i yi J 

calculating the slope and offset of the line of best fit, the m„- — - ----- ■ 

sliding window uses only the PCR sample values within a B " 

window representative of a predetermined time interval. The so 
phrase "sliding window" may be defined as follows. The 

window includes a predetermined number of recently For a window size of n samples, the slope of the line of best 

received sample values and slides in the sense that it is fi t upon analysis of a new PCR sample value can be 

shifted to include each newly-received sample value. calculated using the following recursive equation. 

Sliding, or shifting, may be implemented by adding the 55 

newly-received sample value to the sample values in the (8) 

window and discarding the oldest sample value. Thus, only m "** ""js^T" " 
recently-received PCR sample values arc used in calculating 

the slope and offset. The window terminates with the current , ^ /„ „ \ "v* ^ <v r \"v*„ ^ 

PCR sample value and timestamp. 60 «-* Wfc 

The size of the window may be varied within the scope of (/> + iX*»+*>w-**y0 

the present invention. A preferred window size includes / n+k \ 

about 500 PCR sample values and associated timestamps. + (x* - jc + 0 ( ^"l. x,-(« + + } 

By limiting the number of sample values used in the \ / 

calculation, there is minimal loss in precision due to the 65 

floating point limitation of the computer, and the PCR jitter where n represents the number of sample values in the 

measurement is more accurate. window, k is a positive integer that represents the number of 



n ( n V 
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the window as it slides and i is an index. In Equation (8), 
A n ^ k _ J and B M+Jt _ 3 represent the numerator and denominator, 
respectively, from the previous calculation and in the case 
where k=l, they represent the numerator and denominator, 
respectively, of Equation (7). Similarly, the variance (j) 5 
of the PCR jitter upon analysis of the next sample value, 
given the previously calculated variance Vn (j), is given by 
the following recursive equation. 

v (»-W.tf) <J»*i-j»Y (9) 10 

n « + l 

The PCR jitter variance is calculated in accordance with 
Equation (9) over all sample values, not just the sample 
values within the window. 

A flow chart of a second embodiment of a routine for 15 
determining PCR jitter is shown in FIG. 8. In step 300, a 
PCR sample value is extracted from the transport stream and 
is assigned a timestamp from the local clock. The received 
PCR sample value and the timestamp from the local clock 
are normalized in step 302. Steps 300 and 302 may be 20 
performed as described above in connection with steps 40 
and 42, respectively, of FIG. 2. In step 304, the new 
normalized PCR sample value and timestamp are written in 
a list in memory. The list is an array of sample pairs, 
including the normalized PCR sample values and times- 25 
tamps. A maximum number of sample pairs in the list is 
determined by the window size. In step 306, a determination 
is made whether a sufficient number of PCR sample values 
has been received to determine the line of best fit within a 
given accuracy. When a sufficient number of PCR sample 30 
values has not been received, the process waits for more 
PCR sample values to be received. When a sufficient number 
of PCR sample values has been received, the slope and offset 
of the line of best fit are determined in step 310. The slope 
is determined in accordance with Equation (8), and the offset 35 
is determined in accordance with Equation (4), where the 
summations are performed over the sample values that fall 
within the window. The PCR jitter is determined in step 312 
in accordance with Equation (5). Then the PCR jitter vari- 
ance may optionally be determined in step 314 in accor- 40 
dance with Equation (9). A display of jitter and jitter 
variance may be generated in step 320, as described above 
in connection with step 56 of FIG. 2. In step 324, a 
determination is made as to whether the number sample 
values received is equal to or greater than the number of 45 
sample values in the window. When the number of sample 
values received is less than the window size, the process 
returns to step 300 to wait for another PCR sample value. 
When the number of sample values received is equal to the 
window size, the oldest PCR sample value and timestamp 50 
are removed from the list in step 330. The process then 
returns to step 300 to wait for another PCR sample value. 

While there have been shown and described what are at 
present considered the preferred embodiments of the present 
invention, it will be obvious to those skilled in the art that 55 
various changes and modifications may be made therein 
without departing from the scope of the invention as defined 
by the appended claims. 

What is claimed is: 

1. A method for measuring program clock reference 60 
(PCR) jitter in a data stream containing digitally encoded 
information and PCR samples, comprising the steps of: 
a) providing a test instrument including a line interface 
circuit for receiving the data stream, a digital computer 
connected to the line interface circuit for processing the 65 
data stream and a local clock for generating times- 
tamps; 
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b) said digital computer extracting a PCR sample value 
from the data stream, said PCR sample value repre- 
senting a count of a system time clock used to generate 
said data stream; 

c) said digital computer obtaining a timestamp from said 
local clock and associating the timestamp with said 
PCR sample value, said timestamp representing a time 
of arrival of said PCR sample value; 

d) said digital computer determining a reconstructed 
system time clock from a set of values comprising said 
PCR sample value and said timestamp, and previously- 
received PCR sample values and associated timestamps 
within a window representative of a predetermined 
time interval that terminates with said PCR sample 
value and said timestamp; and 

e) said digital computer determining a PCR jitter value for 
said PCR sample value by determining a deviation of 
said PCR sample value from said reconstructed system 
time clock, said PCR jitter value representing a varia- 
tion in arrival time of said PCR sample value with 
respect to said reconstructed system time clock. 

2. A method as defined in claim 1 wherein the step of 
determining the reconstructed system time clock includes 
determining a line of best fit to said set of values within said 
window, said line of best fit representing said reconstructed 
system time clock. 

3. A method as defined in claim 2 wherein the step of 
determining the line of best fit includes determining a slope 
and an offset of said line of best fit. 

4. A method as defined in claim 1 wherein said data stream 
comprises an MPEG-2 transport stream. 

5. A method as defined in claim 1 wherein the step of 
determining the reconstructed system time clock includes 
performing a recursive linear regression calculation on said 
set of values within said window. 

6. A method as defined in claim 1 further including 
repeating steps b)-e) for a plurality of PCR sample values 
from said data stream to provide a plurality of PCR jitter 
values in real time. 

7. A method as defined in claim 6 wherein the step of 
determining the reconstructed system time clock for each 
new PCR sample value includes shifting said window to 
include said new PCR sample value. 

8. A method as defined in claim 6 further including the 
step of determining the variance of said PCR jitter. 

9. A method as defined in claim 6 further including the 
steps of providing said test instrument with a display unit 
and displaying said plurality of PCR jitter values on said 
display unit. 

10. A method as defined in claim 6 further including the 
steps of providing said test instrument with a display unit 
and displaying said plurality of PCR jitter values in a 
histogram on said display unit, 

11. A test instrument for measuring program clock refer- 
ence (PCR) jitter in a data stream containing digitally 
encoded information and PCR samples, comprising: 

a line interface circuit for receiving the data stream; 

a local clock for generating timestamps; 

a digital computer connected to the line interface circuit 

for processing the data stream, said digital computer 

comprising: 

means for extracting a PCR sample value from the data 
stream, said PCR sample value representing a count of 
a system time clock used to generate said data stream; 

means for obtaining a timestamp from said local clock 
and associating the timestamp with said PCR sample 
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value, said timestamp representing a time of arrival of 
said PCR sample value; 

means for determining a reconstructed system time clock 
from a set of values comprising said PCR sample value 
and said timestamp, and previously-received PCR 
sample values and associated timeslamps within a 
window representative of a predetermined time interval 
that terminates with said PCR sample value and said 
timestamp; and 

means for determining a PCR jitter value for said PCR 
sample value by determining a deviation of said PCR 
sample value from said reconstructed system lime 
clock. 

12. A test instrument as defined in claim 11 wherein said 
means for determining the reconstructed system time clock 
includes means for determining a line of best fit to said set 
of values within said window, said line of best fit represent- 
ing said reconstructed system time clock. 

13. A test instrument as defined in claim 12 wherein said 
means for determining the line of best fit includes means for 
determining a slope and an offset of said line of best fit. 

14. A test instrument as defined in claim 11 wherein said 
data stream comprises an MPEG-2 transport stream. 

15. A test instrument as defined in claim 11 wherein said 
means for determining the reconstructed system time clock 
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includes means for performing a recursive linear regression 
calculation on said set of values within said window. 

16. A test instrument as defined in claim 11 wherein said 
digital computer includes means for updating said recon- 

5 structed system time clock and means for determining a 
PCR jitter value for each of a plurality of new PCR sample 
values received in said data stream to provide a plurality of 
PCR jitter values. 

17. A test instrument as defined in claim 16 wherein said 
10 means for updating said reconstructed system time clock 

includes means for shifting said window to include each of 
said new PCR sample values. 

18. A test instrument as defined in claim 16 wherein said 
digital computer further includes means for determining the 

15 variance of said PCR jitter. 

19. A test instrument as defined in claim 16 wherein said 
digital computer further includes a display unit and means 
for displaying said plurality of PCR jitter values on said 
display unit. 

20 20. A test instrument as defined in claim 16 wherein said 
digital computer further includes a display unit and means 
for displaying said plurality of PCR jitter values in a 
histogram on said display unit. 

* * * * * 
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